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ABSTRACT 

Context. In several instances chemical abundances of main-sequence and giant stars are used simultaneously under the assumption 
that they share the same abundance scale. This assumption, if wrong, might have important implications in different astrophysical 
contexts. 

Aims. It is therefore crucial to understand whether the metallicity or abundance differences among dwarfs and giants are real or are 
produced by systematic errors in the analysis. We aim to ascertain a methodology capable of producing a consistent metallicity scale 
for giants and dwarfs. 

Methods. To achieve that, we analyzed giants and dwarfs in the Hyades open cluster, under the assumption that they share the same 
chemical composition. All the stars in this cluster have archival high-resolution spectroscopic data obtained with HARPS and UVES. 
In addition, the giants have interferometric measurements of the angular diameters. We analyzed the sample with two methods. The 
first method constrains the atmospheric parameters independently from spectroscopic method. For that we present a novel calibration 
of microturbulence based on 3D model atmospheres. The second method is the classical spectroscopic analysis based on Fe lines. We 
also tested two different line lists in an attempt to minimize possible non-LTE effects and to optimize the treatment of the giants. 
Results. We show that it is possible to obtain a consistent metallicity scale between dwarfs and giants. The preferred method should 
constrain the three parameters T^-g, log g, and ^ independent of spectroscopy. A careful selection of Fe lines is also important. In 
particular, the lines should not be chosen based on the Sun or other dwarfs, but specifically to be free of blends in the spectra of 
giants. When attention is paid to the line list, the classical spectroscopic method can also produce consistent results. In our test, the 
metallicities derived with the well-constrained set of stellar parameters are consistent independent of the line list used. Therefore, for 
this cluster we favor the metallicity of +0.18±0.03 dex obtained with this method. The classical spectroscopic analysis, using the line 
list optimized for the giants, provides a metallicity of •f0.14±0.03 dex, in agreement with previous works. 

Key words. Stars: abundances - Stars: fundamental parameters - Stars: late-type - Techniques: spectroscopic - Galaxy: evolution 
- Planets and satellites: formation 


1. Introduction 


The determination of chemical abundances in stars by spectral 
synthesis or curve of growth is a rather well-established tech¬ 
nique, largely available to most astronomers. A more detailed 
view of stellar abundances has become possible since the advent 
of 8-10m class telescopes coupled with high-efficiency spectro¬ 
graphs, which produced high quality spectra for many stars. As 
a consequence, uncertainties in abundance analyses are mainly 
dominated by systematic errors rather than by observational lim¬ 
itations regarding spectral resolution and/or signal-to-noise ratio 
(S/N). In this context, we should recall that the determination 


* Based on data obtained from the ESO Science Archive Facility. 
The observations were made with ESO Telescopes at the La Silla and 
Paranal Observatories under programmes ID 070.D-0356, 088.C-0513 
and 070.D-042L 

** e-mail: leticia@dfte.ufrn.br 


of stellar abundances requires a number of assumptions, some 
of which might not be valid for stars in different evolutionary 
stages, such as dwarfs and giants. 

Although cool dwarfs and giants are often analyzed follow¬ 
ing a single methodology, the differences in their photospheric 
properties might introduce distinct systematic effects on the fi¬ 
nal atmospheric parameters derived during the analysis. An ex¬ 
ample is the large experiment with multiple analysis pipelines 
conducted within the Gaia-ESO Spectroscopic Survey. In this 
survey, it was clearly demonstrated that the multiple analysis 
methodologies per form differently in dist inct regions of the pa¬ 
rameter space (see ISmilianic et £0120141) . It is thus important 
to understand if differences on metallicities and abundances be¬ 
tween giants and dwarfs are real or produced by systematic er¬ 
rors in the analysis. The interpretation of some open questions 
in astronomy depends on these kinds of comparisons. 
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For instance, metallicity plays a role in one of the most ac¬ 
cepted scenarios in planet for mation theories, which is the core 
accretion scenario, (see, e.g.. iPollack et al.l[T99^ . Concern¬ 
ing the gas-giant planet formation, many works have shown 
that main-sequence FGK-type stars hosting giant planets (> 1.0 
Mj) are usually metal-rich w hen compared with their counter- 
parts_without_gian^^lanets(|Fischer_^^^lentilj200^Gonza- 
lez Il997l : ISantos et al.ll2004l) . However, a giant planet versus 
metallicity correlation is not well established for evolved stars. 
iPasQuini et al.l (l2007h showed that giant stars with planets do 
not tend to have high metallicities. A similar result w as found 
by othgrworkgj|Gh^ezzieLa]J|2010LlMortie reLa]J|2013LSchuler 


et al. l2() 0.5|: [Takeda et alj ^OOSh . but not bv iHekker & Melendea 

(l2007h : lReffert et al.l (1201 5h \^o also argued in favor of a planet 
metallicity correlation in giants. 

The situation for giants is more complex because of many 
factors as, for example, the limited sample size of giants host¬ 
ing planets, inhomogeneities in the planetary properties, and the 
higher masses of giant stars in comparison to dwarfs. Indeed, 
there is a stellar mass vs. planet vs. metallicity correlation when 
we consider the giant planet occurrence frequency in evolved 
stars. Some authors have reported that there is a trend toward 
highemnetalIicities^orstarswithrnasses_aboveJ_j5MQ_^Johnson 
et al. l2010l : lMaldonado et alj|2013h . Understanding the metallic - 
ity distribution of stars with planets would provide an important 
constraint for planet formation theories. 

The comparison between the metallicity scale of giants and 
dwarfs is also important in studies of Galactic chemical evo¬ 
lution. For instance, the metallicity distribution (MD) of the 
Galactic bulge provides clues about how similar bulge stars are 
to thin and thick disk stars. As a consequence, it aids the deter¬ 
mination of the bulge ag e. Initial studies of the bulge were fo¬ 
cused on giant stars (e.g.. lMc William & Richll99^IZoccali et al.l 
I 2 OO 6 II because these are intrinsically brighter objects. Later, the 
microlensing technique was used to observe dwarfs and sub¬ 
giants in the Galactic bulge. This reveal ed discrepancies be¬ 
tween the MD of bulge g iants and dwarfs dBensbv et al.ll20lll 
l2010l : l^hen et all 120081) More recent stu dies, however, show 
a better agreement between the two MDs (I Anders et alj|201^ : 
iBensbv et al.l 1201^ . although the authors recognize some ev¬ 
idence of a bias toward the high-metallicity tail of the giants’ 
distribution (iHill et alj|201 IfciTavlor & Croxallll2005h . Indeed, it 
is difficult to analyze giants in the high-metallicity regime since 
their spectra are severely affected by blending and molecules fea¬ 
tures because of their cool atmospheres. A better picture of the 
MD of giants requires a full knowledge of the source of this bias 
in the metal-rich regime. 

Methodological limitations may affect the analysis of giants 
and dwarfs differently, e.g., how realistic are the adopted at¬ 
mospheric models, continuum normalization, and atomic and 
molecular constants used. The use of the same line list can be 
a challenge since the intensity of the spectral lines is different 
in these objects. Moreover, departures from the local thermo¬ 
dynamical equilibrium (LTE) are particularly im portant for low- 
gravity and low-metallicity stars (lAsDlundl200.5h . Spectroscopic 
gravities derived by the ionization equilibrium may be unreli¬ 
able because of Fei overionization. High-temperature dwarfs 
(Teff > 6000 K) ma y also suffer from non-LTE effects ( Mashon- 
kina et al. l2010l) . Additionally, the gf values need to be very 
accurate for giants since a differential analysis with respect to 
the Sun does not cancel out uncertainties in these constants. 

Open clusters are the optimal sites to evaluate the limitations 
of abundance analyses when applied to giants and dwarfs. They 
are made of stars with basically the same distance, and in gen¬ 


eral, it is reasonable to assume that all stars share the same over¬ 
all chemical composition, except for the element s affected by 
mixing in giants (see, e.g ., ISmilianic et al.l 120091 : ISmith et al.l 
I 2 OI 3 I : l^keda et al]l2008ll . Stars in an open cluster also share 
the same age, outlining a common isochrone curve in the HR 
diagram, making it easier to constrain stellar parameters in these 
environments. 

Eew studies so far have attempted the simultaneous analysis 
of giants and dwarfs in open clusters, aiming to explore possi¬ 
ble discrepancies in the abundance pattern s between these two 
classes of objects. iPasauini et al.l (l2004l) investigated dwarfs 
and giants of the intermediate-age cluster IC 4651, and found, 
in general, exc ellent agreement b etween stars of different evolu¬ 
tionary status. iPace et akl (1201 Ol) investigated the metallicity of 
five dwarfs and three giants in two open clusters and found differ¬ 
ences in the metallicity of up to 0.10 dex for one of the clusters. 
They also pointed out enhancements of sodium, aluminium, and 
silicon for^he_giantSj^A_simite_snjd^_^as_£erformed_^San- 
tos et al. (I2012il2009ll . who investigated the abundance pattern 
of several open clusters performing a simultaneous and homo¬ 
geneous spectroscopic analysis. These authors noticed that the 
discrepancy on the metallicity scale of giants and dwarfs belong¬ 
ing to the same cluster may depend on the line list used. The 
explanation about the source of these differences is still a matter 
of investigation. 

In this work, we chose the Hyades open cluster to perform a 
simultaneous and homogeneous study of giant and dwarf stars. 
Our aim is to define an analysis method that can deliver metal¬ 
licities in a consistent scale for both types of stars. Once this 
method is tested and established, it can be applied to different as- 
trophysical problems, such as the comparison of MDs of planet 
host dwarfs and giants. We choose the Hyades, as atmospheric 
parameters can be constrained by other methods than the classi¬ 
cal spectroscopic analysis (see below), and because this cluster 
is young enough that it should be free of atomic diffusion effects. 
The paper is organized as follows. Section |2] reviews the main 
properties of the cluster, while Section [3]presents a description 
of the data. Section|4]is dedicated to the analysis and Section|5] 
to the discussion of the results. In Section |6] we draw our final 
conclusions. 


2. Hyades: The benchmark test 

The Hyades is a re latively young cluster , with an estimated age 
of ~ 625 + 50 Myr (iPerrvma n et al.ll 9^^. and is the closest open 


cluster to the Sun (~ 46.5 pc, van Leeuwenl200^ . Spectroscopic 


studies of EGK-t ype dwarfs in the Hyades find a metallicity of 
about -1-0.13 dex ( Boesgaa rd & Eri ^ j_990HCavrel et "^119851 : 
iPaulson et al.ll2003 : Schuler et al. 2006ll . Regarding the giants, 
the metallicity values range from -I-O.IO up to -(-0.20 dex, where 
this scatter is usually attributed to the star HIP 20455, a spec tro- 
scopic binary (ICarrera & Pancin3201 lUSchuler et al.ll20()6ll . 

The slightly over-solar metallicity provides a safe regime to 
test the classical spectroscopic analysis. At this regime, depar¬ 
tures from LTE are not expected to be significant for Ee i. Also, 
close to the solar metallicity, the mean temperature stratifica¬ 
tion is close to the radiative equilibrium expectation, and there¬ 
fore, the difference between I D and 3D mod el atmospheres is 
expected to be relatively small (lAsDlundll2005h . 

A relevant advantage of the stars in the Hyades is that their 
atmospheric parameters can be very well constrained. Erom the 
standpoint of the giants, interferometric measurements of the an¬ 
gular diameters are available, which enab les the direct determi - 
nation of absolute effective temperatures (iBovaiian et al.ll2()()^ . 
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Furthermore, precise Hipparcos parallaxes dvan Leeuwenll200^ 
are available. Thus, a reliable determination of the surface grav¬ 
ities is possible. Finally, a large amount of spectroscopic data 
is available for the numerous dwarfs and all four giants of this 
cluster. All these make the Hyades the optimal benchmark to test 
the limitations of the classical spectroscopic analysis method. 

The Hyades is a lso a target of planet searches using the ra dial 
velocity technique (ICochran et al.i20()^lPaulson et al 1200411 . So 
far, one of the four giants in the Hyades was reported to host a 
giant planet. The clump giant HIP 20889 has a l ong period planet 
(~ 594 days) with ~ 7.6 Mj (ISato et al.ll2007h . More recently, 
lOuinn etdl (1201 4l) reported the discovery of the first hot Jupiter 
orbiting a K dwarf HD 285507 in this cluster. 

Our sample was selected as f ollows. The dwarf s were se¬ 
lected primarily from the list of Ivan BuerenI (Il952h and then 
cross-ch ecked with the r e liable sample of cluster members de¬ 
fined by iPerrvman et al.l ( 1998h. Three addition al dwarfs were 
selected exclusively from iPerrvman et ^ (Il998h to complement 
the cooler end of our sample (these are the stars without the vB 
number in Table [T]i. Among them, one cool dwarf (HIP 13976) 
appears slightly away from the cluster’s main sequence. How¬ 
ever, t his object was classified as a reliable member of the clus¬ 
ter by IPerrvman et al.l (Il998h and has radial velocity and dis¬ 
tance fully compatible with the cluster distribution. Moreover, 
this starjs_presentjn^tun^_ana]vses_ofdiecluster_(e;^Paulson 
et al. 120031: LSchiiler et al 1120061 : 1 Yon g et alJl2nn4h . 

SeverdstudiesjnvestigatedbinariesjnJheJJ^ades_{Stefanik 
& Latham n992LI 19851) and our stars were chosen to avoid binary 
systems. Among the giants, we excluded HIP 208 85, as it is a 
spectr oscopic binary (SBl) with a blue companion. iTorres et al.l 
(Il997h estimated that the secondary star contributes about 3% of 
the flux of the primary. An accurate abundance analysis should 
take this contribution into account. Finally, whenever possible, 
we chose stars that have been s tudied in previous works for com¬ 
pariso n purposes (in particular. lPaulson et al.l20()^ : ISchuler et al.l 
l2006l) . Table [T] presents the basic data of the 14 dwarfs and the 
three giants selected for our sample. Their position in the CMD 
is shown in Fi g.[T| In this figure, the magnitude V and the (B - V) 
color are from iPerrvman et aH(ll998l) . We highlight that our sam¬ 
ple encompasses stars in a large range of effective temperatures 
(4700 K < Teff < 6200 K). This facilitates the investigation of 
possible systematic effects from the analysis as a function of this 
parameter. 


3. Observations and data reduction 

All spectra used in this work were downloaded from the Euro¬ 
pean Southern Observatory (ESO) science archive facilitjQ. The 
giants were observed during ESO program 088.C-0513 with the 
HARP^ffligh^ccuragy^gdialvelocit^^lanet^eMcher^Mayor 
et al. l2003h high-resolution spectrograph IR = 110 000). situated 
at the 3.6 m telescope in La Silla, Chile. The spectral range 
covers from 3800 to 6900 A, with a small gap between 5300- 
5330 A because of the arrangement of the CCD mosaic. The 
average S/N @6109 A is about 400. All the spectra were repro¬ 
cessed by the last version of the HARPS pipeline (Data Reduc¬ 
tion Software version 3.5). We only carried out the Doppler cor¬ 
rection and the continuum normalization with standard IRAfQ 
routines. 


* http://archive.eso.org/eso/eso_archive_adp.html 
^ Image Reduction and Analysis Facility - IRAF provided by Associa¬ 
tion of Universities for Research in Astronomy- AURA, EUA. 
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Fig. 1. Color magnitude diagram (CMD) of the Hyades. The large 
black triangles correspond to the sample selected for our analysis. The 
small gray dots are the oth er cluster members. T he magnitude V and 
the (B — V) color are from IPerrvman et all d 19981) and all stars shown 
fulfill the membership criterion of this work. 


The dwarfs were observed during ESO program 70D- 
0356 with UVES (Ul traviolet and Visual Echelle Spectrograph, 
iDekker et al.l l200fll) at the 8.2m Kueyen telescope of the 
VET (Very Large Telescope) with and spectral resolution of 
R ~ 60 000. Eor consistency, we used spectra acquired with the 
same instrument configuration for all the observations: slit width 
of 0.8" and central wavelength at 580 nm in the red arm. The 
wavelength coverage is from 4780-6800 Awith a gap between 
5750-5830 A because of thee arrangement of the CCD mosaic. 
The S/N varies from 150 to 400. Da ta reduction was carrie d out 
with the ESO Reflex environment (lEreudling et al.ll2013l) ver¬ 
sion 2.3, with the exception of Doppler correction, which was 
performed with IRAE routines. 

We adopt two spectra as solar proxies. Eor comparison with 
the giants, we use the reflected spectrum of the Jupiter’s moon 
Ganymed^l, which was obtained with the same HARPS config¬ 
uration as the spectra of the giants. Eor the comparison with the 
dwarfs, we use the UVES solar spectruirQ , which was obtained 
with the moonlight illuminating the slit. Both spectra have, on 
average, a S/N > 300. 


4. Analysis 

The atmospheric parameters (effective temperature, Tgff, surface 
gravity, log g, and microturbulence, ^) of the sample stars were 
determined using two different methods. In this way, we can 
compare the final metallicity scale obtained under different as¬ 
sumptions. The first method, hereafter Ml and described in Sec- 
tion l4.3l takes constraints into account that do not depend on the 
classical spectroscopic method, but that are still fine-tuned us¬ 
ing Eei and Eeii lines. The second method, hereafter M2 and 

^ http://WWW.eso.org/sci/facilities/lasilla/ 
instruments/harps/inst/monitoring/sun.html 

http: //WWW. eso. org/observing/dfo/quality/UVES/ 
pipeline/solar_spectrum.html 
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Table 1. Selected_obiectsJnJhe_^l^ades_clusterj_TheJettersJnJhe_columns_^efe£^o;_^a]_^an^uerenjiumber^_^b]_^jnagnitude_^o^ & 

Knuckles (ll9^, exce pt for the stars with number marks; (c) Hipparcos parallax and its respective standard error (mas); (d) Radial velocity from 
iPerrvman et al.l ( ll998h and its respective error; (e) adopted mass in solar units (see Section lTSl for details). 


Giants 


HIP 

vB 

Spec. Type 

V 

RA(J2000) 

Dec(J2000) 

(?r ± CTn ) 

(K ± o-vf 

Mo 

S/N 







(mas) 

(km s“‘) 


@ 609 nm 


(a) 


(b) 



(c) 

(d) 

(e) 


20205 

28 

KOIII 

3.66 

04:19:47.6 

-h15:37:39.5 

21.17+1.17 

+39.28+0.11 

2.48 

400 

20455 

41 

KOIV 

3.77 

04:22:56.1 

+ 17:32:33.0 

21.29+0.93 

+39.65+0.08 

2.48 

380 

20889 

70 

KOIII 

3.52 

04:28:36.9 

+ 19:10:49.5 

21.04+0.82 

+39.37+0.06 

2.48 

440 

Dwarfs 

13976 


K2.5V 

‘7.95 

03:00:02.8 

+07:44:59.1 

42.66+1.22 

+28.35+0.18 

0.83 

220 

16529 

4 

G5D 

8.88 

03:32:50.1 

+23:41:31.9 

22.78+1.26 

+32.72+0.17 

0.87 

210 

18946 


K5D 

M0.13 

04:03:39.0 

+ 19:27:18.0 

23.07+2.12 

+36.93+0.26 

0.75 

150 

19098 


K2D 

2 9.29 

04:05:39.7 

+ 17:56:15.7 

19.81+1.39 

+37.61+0.05 

0.88 

160 

19148 

10 

GOV 

7.85 

04:06:16.1 

+ 15:41:53.2 

21.41 + 1.47 

+38.04+0.17 

1.08 

310 

19781 

17 

G5V 

8.46 

04:14:25.6 

+ 14:37:30.1 

21.91 + 1.27 

+39.24+0.06 

0.97 

290 

19793 

15 

G3V 

8.09 

04:14:32.3 

+23:34:29.8 

21.69+1.14 

+38.21+0.23 

1.01 

320 

19934 

21 

G5D 

9.15 

04:16:33.5 

+21:54:26.9 

19.48+1.17 

+38.46+0.19 

0.92 

230 

20130 

26 

G9V 

8.63 

04:18:57.9 

+ 19:54:24.1 

23.53+1.25 

+39.58+0.06 

0.93 

330 

20146 

27 

G8V 

8.46 

04:19:08.0 

+ 17:31:29.1 

21.24+1.32 

+38.80+0.08 

0.94 

300 

20899 

73 

G2V 

7.85 

04:28:48.3 

+ 17:17:07.7 

21.09+1.08 

+39.37+0.06 

1.06 

470 

21112 

88 

F9V 

7.78 

04:31:29.3 

+ 13:54:12.5 

19.46+1.02 

+40.98+0.31 

1.13 

320 

22422 

118 

F8D 

7.74 

04:49:32.1 

+ 15:53:19.5 

19.68+0.96 

+42.04+0.14 

1.10 

320 

22566 

143 

F8D 

7.90 

04:51:23.2 

+ 15:26:00.5 

17.14+1.00 

+42.92+0.19 

1.17 

250 


Notes. Additional sources of V magnitudesdli lKoen et al.U201Cl) . t2~) I.Tohnson et al.l(ll962h . 


described in Section 14.41 is the classical spectroscopic analysis 
where the parameters are determined using the Fe i and Fe it ion¬ 
ization and excitation equilibria. In addition, in the implementa¬ 
tion of each method, we make use of two different line lists. We, 
therefore, derived four sets of atmospheric parameters for each 
star. Before describing the two methods, we present the two line 
lists and discuss the measurement of equivalent widths (EWs). 

4.1. Line lists selection 

The two line lists that we adopted were assembled with two 
goals. First, we aim to minimize non-LTE effects as it will 
differentially affect stars of different gravities and temperatures. 
Therefore, the first list contain a set of lines for which non-LTE 
effects would be minimiz ed or neglected, at least in the metallic- 
ity regime of the Hvades. flVlashonkina et al.l(l201 ih evaluated the 
non-LTE line formation of the two ions of iron in cool reference 
stars, some of them, with metallicity comparable to the Hyades. 
In that study, the authors concluded that non-LTE effects are vir¬ 
tually negligible for Eeii lines (i.e., they affect the abundances 
by less than 0.01 dex) and are very small for Ee i lines, for stars 
of metallicity slightly higher than the Sun. We selected a to¬ 
tal of 42 Eei and 15 Eeii, from their line list, among weak to 
moderately strong transitions, which were well isolated and as 
free as possible from blending features. We excluded 14 lines of 
the original list for which the measured EW in the Sun was in 
the sat urated regime of the curve of growth. iMashonkina et al.l 
(1201 ih claimed that the accuracy of their iron abundances might 
be affected by the uncertainties in the gf values that were used. 
Their gf values were obtained from experimental measurements 
collected from different papers. We improved some of the gf 


values in this line list with more recent determinations that were 
kindly provided by Maria Bergemann (private communication). 
Table lA.fl lists the selected lines, their atomic data and, in addi¬ 
tion, the equivalent widths and individual abundances obtained 
for the solar reflected spectrum of Ganymede using ML This 
line list is hereafter referred to as MASH. 

The second list was chosen to be suitable for the analysis 
of giants. Often, line lists assembled for the analysis of the Sun 
may not be optimized for giants because of the more pronounced 
spectral transitions present in these stars. We have used a line 
list with transitions carefully chosen to avoid blends in giants 
and with accurate gf values determination provided through the 
courtesy of Dr. Martin Asplund (private communication). This 
list includes a total of 34 Fe i and 7 Fe ii transitions and here¬ 
after is referred to as ASPL. Table IA.2I lists the lines, as well 
their atomic data, and the equivalent widths and individual abun¬ 
dances for the solar reflected spectrum of Ganymede using ML 
There are 15 lines in common between the MASH and ASPL 
lists. Nevertheless, the abundances in the Sun are sometimes 
distinct because of the different atomic data adopted in each list. 
This difference is greater than 0.05 dex for about 20% of the 
lines, while for the remaining it is about ~ 0.02 dex on average. 

4.2. Equivalent width measurements 

We used the code ARES (Automatic Routi ne for line Equivalent 
widths in stellar Spectra JSousa et al.l2007h to perform automatic 
measurements of the EWs of the Fei and Eeii lines. This code 
applies a Gaussian fit to the profile of the absorption lines, tak¬ 
ing a local continuum into account, which is determined on the 
basis of a set of input parameters provided by the user. We tested 


Article number, page 4 ofll9l 

























Dutra-Ferreira, L. et al.: Consistent metallicity scale for cool dwarfs and giants 


different combinations of input parameters for both UVES and 
HARPS spectra, and chose those which visually produced the 
best fits to the line profile. The best fits were determined by a 
visual check and by comparing the EWs computed with ARES 
and the EWs computed with IRAE for a given set of stars. These 
were the same stars selected for the EWs comparison between 
IRAE and ARES (see below). With this approach, the input fit¬ 
ting parameters for ARES were optimized. The spectral reso¬ 
lution we use is such that, in general, the instrumental profile 
dominates the observed profile, and therefore, a Gaussian fit can 
well reproduce the observed profile of the lines. Additionally, 
we removed all EWs < 5 mA from our analysis to avoid lines 
severely affected by noise or uncertainties related to the contin¬ 
uum fit, therefore, with larger relative errors in the EW measure¬ 
ment. We also removed EWs >120 mA, to avoid the flat part of 
the curve of growth, in which a Gaussian fit may not adequately 
reproduce the observed line profiles. 

The quality of the automatic measurements is comparable to 
the manual method obtained with the task splot of IRAE for the 
majorit^_ofthe^^^^S_an^^VES_S£ectru_as_shown_^^Sousa 
et al. (l2007h . Notwithstanding, we repeated this comparison 
with the line lists used in our work to verify if the manual EWs 
obtained with splot are comparable with those computed using 
ARES for both giants and dwarfs. 

Eigure |2] shows the results of such comparison. In ad¬ 
dition to the solar spectrum, we select as examples the hot 
dwarf HIP 22566, where the continuum placement is more 
challenging because of the enhanced rotation, and the giant 
HIP 20205, where we have fewer good lines available for the 
analysis. There is in general a good agreement between the au¬ 
tomatic and the manual measurements. The mean difference be¬ 
tween the measurements is 0.717 + 1.505 mA for Ganymede, 
1.907 + 1.933 mA for HIP 20205, and 1.074 + 2.413 mA for 
HIP 22566. As seen in Pig. the slope and the intercept values 
of the linear regressions between the two sets of EWs are close 
to one and zero, respectively. 

Although the difference between the two sets of EWs is 
small, it is not negligible, but the comparison indicates the same 
trend for all stars tested, i.e, the EWs obtained with splot are sys¬ 
tematically higher than those obtained with ARES. To quantify 
the effect of the EW measurement differences on the final metal- 
licities, we computed the abundances for the stars shown in Pig|2] 
using both the automatic and the manual measurements. The 
mean difference of the metallicities is ~ 0.03 dex for the Sun, 
and ~ 0.04 dex for both HIP 22205 and HIP 22566. Based on 
the small differential variation found in the derived abundances, 
we therefore decided to measure all the EWs using ARES. We 
visually inspected the fits for each line, and removed from our 
analysis those fits that were judged of poor quality. Bad line fits 
were mainly attributed to poor estimates of the local continuum 
regions or lines that were significantly affected by noise. 

4.3. Method M1: Well-constrained parameters independent 
of spectroscopy 

The first method used to constrain the atmospheric parameters 
makes use of input values that have been determined indepen¬ 
dent of the classical spectroscopic method. With this new ap¬ 
proach, we expect to be able to investigate which are the pos¬ 
sible systematics that can arise from the classical spectroscopic 
analysis. 

To r the red giants, we ad opted the direct determinations of 
Teff bv iBovaiian et al.l (l2009ll . Angular diameters were obtained 


with long-baseline optical interferometry and transformed into 
linear radii using the Hipparcos parallaxes. The direct estima¬ 
tions of Teff are calculated combining these radi i and the bolo- 
metric flux of the star (see iBoehm-Vitenselll 98 iL for the defini¬ 
tion of direct measurements of for EGK stars). The bolomet- 
ric fluxes were determined using t he bolometric corrections from 
lAllende Prieto & LambertI (Il999ll and assuming Mbol.o = 4.74. 
Effective temperatures derived by this method can reach an ac¬ 
curacy of 1% and are the state of the art in Teff determinations. 

Eor the dwarfs, as interferometric angular diameters are not 
available, we obtained from photometric calibrations derived 
with the InfraRed Elux Method (IREM). Although the IREM 
provides semidirect estimations of the temperature scale is 
almost m odel independent. We use the JHKs magnitudes from 
2MASS (ISkrutskie et al.l 2 006ll an d the {B - V) color from the 
Hipparcos catalog (I van LeeuwenI 200% . We adopt the calibra¬ 
tions from ICasagrande et al.l (1201 Oh to derive photometric tem¬ 
peratures using the colors (V-J), (V -H), (V-Ks), and (B-V). 
The adopted photometric Teff was calculated with an average 
of the four different temperatures estimates weighted by the er¬ 
rors of the calibrations. Effective temperatur es derived by this 
metho d reach and accuracy of a few percent (ICasagrande et al.l 

IMH- 

Since the Egfr of the giants were derived from interferomet¬ 
ric measurements and those of the dwarfs were derived from 
IREM, it is interesting to evaluate how well the two methods 
agree fo r the giants. The IREM T^g for the giants are shown in 
Table |2] Ida Silva et al.l (l2006ah showed that 2MASS colors are 
unsuitable to determine Egfr of bright stars. Eor this reason, we 
used only the (B-V) color combined with the calibrations of 
lG^,:^lez Hernandez & B onifaciol(l2009ll . l^mfrez & Melendej 
(l2005h arid Alonso et al. (1 19961) . vvhich are more appropriated for 
evolved stars. Interferometric effecti ve temperatures (7’eff(int)) 
obtained from iBovaiian et al.l (l2009l) are also shown for com¬ 
parison in Table |2] The average IREM T^g are in a excellent 
agreement with the interferometric values within less than 50 K, 
as shown in the last column of Table |2] A difference between 
individual IREM and interferometric Tf.g as large as 100 K is 
fo und only for the star HIP 20889 when u sing the calibration 
of iGonzalez Hernandez & Bonifaciol (l2009h . which gives values 
systematically lower than the others. 

Surface gravities were determined according to the following 
equation: 

log(g./go) = log(M./Mo)H-41og(TefF./Teffo)-logg(L./Lo). (1) 

Eor the dwarfs, the m asses were computed with the theoretical 
evolutionary tracks of iGirardi et al.l (l2000l) and a Bayesian es¬ 
timation method, which takes the err or of each quantity for the 
mass determination into account (see Ida Silva et al.ll2006bl for 
details). Luminosities were taken from ide Bruiine et al.l (1200 Ih 
and the T^g are those computed as mentioned above. Eor the gi¬ 
ants we a dopted a value of ma ss of 2.48 Mq; the mass of a clump 
giant in a lGirardi e t al.l (1200^ isochrone of 625 Myr and [Ee/H] 
= -1-0.13 (see lSmilianicll201Z for details). 

One free parameter that can not be constrained by our obser¬ 
vational knowledge of the cluster is the microturbulence velocity 
(f). Eor Ml, we fixed ^ adopting predictions obtained with 3D 
atmospheric models. The details on how these values were ob¬ 
tained are described in Section l4.3.1l 

To calculate the metallicities, we used ID-LTE plane-parallel 
atmospheric m odels computed using the Linux vers ion of the 
ATLAS9 code (ISbordonel2005l :ISbordone_e^l.l l2004l) originally 
developed by Kurucz (see, e.g.. lKuru czlll993l ) and adop t ing th e 
ODENEW opacity distribution from Castelli & Kurucd (l2003l) . 
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Fig. 2. Comparison of EWs obtained with ARES and splot. Left: the solar reflected spectrum of Ganymede; middle: the star HIP 20205; and 
right: the star HIP 22566. The one-to-one relation is shown with the dashed line, and the red line represents the linear regression fit in each plot. 
The slope, the intercept, and the rms of the fits are also given in each plot. 

Table 2. IREM obtained with three different calibrations for the giants of the Hyades. TefTfint) obtained from iBovaiian et ^ ( l2009ll are also 
shown for comparison. The < Teff > corresponds to the average of the three IREM T^g. The last column corresponds to the difference between the 
< Teff > and TeffCint). 


HIP 

7'eff(GH09) 

7’eff(AL96) 

reff(RM05) 

< Teff > 

Teff(int) 

ATeff 

20205 

4782 

4874 

4892 

4849 

4844 

5 

20455 

4778 

4870 

4888 

4845 

4826 

19 

20889 

4720 

4810 

4820 

4783 

4827 

-44 


Notes . GH09, AL96 ,and RM05 stands for the calibrations of iGonzalez Hernandez & Bonifacid ll2009h . lAlonso et all (Il996h and I Alonso et al.l 
( I1996I1 . respectively. 


The mixing length parameter adopted was 1.25 and no over¬ 
shooting was considered for both giants and dwarfs. Chemical 
abu ndances of Fe i and Feii were derived using WIDTH pack¬ 
age |^ru^([l993|), under some minor optimizations to facilitate 
handling the input data. 

Our main innovation in method Ml is that the spectroscopy 
independent parameters described above were used as input val¬ 
ues and subjected to a further controlled fine-tuning. The best set 
of parameters for each star was determined as follows. First, we 
allowed the stellar parameters to vary within a range of conser¬ 
vative errors to find the best agreement between the abundances 
of the Fei and Feii lines. The ranges were +50 K for Teff, with 
steps of 10 K; +0.10 dex for log g, with steps of 0.05 dex; and 
+0.10 km s ' for with steps of 0.05 km s *. These ranges were 
chosen because they represent typical error values of the classi¬ 
cal abundance analysis. We tested variations using smaller steps, 
but no major improvement on the final solution was found. Next, 
after reaching a solution, we applied a 2cr-clipping to remove 
lines classified as outliers with respect to the average abundances 
of the previous iteration. The final set of atmospheric parame¬ 
ters is given by looking for the best possible agreement between 
the abundances of Fe i and Fe ii in the 3D-plane - log g - 
The stellar parameters and metallicities derived with this ap¬ 
proach, for giants and dwarfs, with both MASH and ASPL lists, 
are shown in Table [3] 


The main aspect of Ml is that the final parameters shown in 
Table[3]can vary within a very narrow range, constrained by in¬ 
dependent methods. Thus, the errors for the stellar parameters 
in Ml are the uncertainties in the input parameters. An error of 
about ~50 K corresponds to the error of the IRFM calibrations 
and is comp arable with the inter ferometric errors presented in 
the work of (iBovaiian et al.ll200^ see their Table 4). For log g, 
we adopted an error of 0.10 dex because we have very small er¬ 
rors in the parallaxes, thus the main source of error in the gravi¬ 
ties comes from our mass determination. An error of 20% in the 
masses changes the surface gravities by ~ 0.10 dex. Since the 
evolutionary status of the Hyades is well known, we consider 
the gravity error of 0.1 dex as conservative and most likely the 
highest source of uncertainty is in the mass-loss estimate for the 
giants. For ^ we adopted an error of 0.10 km s“*, which corre¬ 
sponds to the uncertainty estimated with the 3D microturbulence 
calibration (see Section 14.3.11) . The dispersion for the Fei and 
Fe II abundances shown in Table |3] correspond to the standard 
deviation of each abundance distribution. 

We also evaluate how the uncertainty in the physical param¬ 
eters of the stars affect the retrieved abundances when using Ml. 
To this end, we assumed the parameters presented in Table [3]for 
the giant HIP 20205 and the dwarf HIP 19148, and then vary 
Teff, log g and f separately for both ASPL and MASH lists. The 
results are given in Table |4] We emphasize that the variations 
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Table 3. Stellar parameters and iron abundances for giants and dwarfs derived according to method Ml. The left side of the table shows the results obtained with the MASH line list and the right 
side the results obtained with the ASPL line list. 






MASH 






ASPL 



HIP 

Teff 

log 9 


logfFcI 

N(Fel) 

logfFcII 

N(Fell) 

Teff 

log 9 


logfFcI 

N(Fel) 

logepell 

N(Fell) 

20205 

4874 

2.61 

1.30 

7.62±0.10 

22 

7.62±0.09 

11 

4884 

2.61 

1.30 

7.59±0.06 

20 

7.59±0.06 

7 

20455 

4876 

2.59 

1.30 

7.59±0.10 

22 

7.62±0.09 

11 

4876 

2.59 

1.35 

7.56±0.07 

20 

7.56±0.06 

7 

20889 

4817 

2.65 

1.35 

7.70±0.09 

18 

7.70±0.10 

10 

4827 

2.60 

1.35 

7.63±0.08 

18 

7.63±0.06 

6 

average 




7.64±0.06 


7.65 ±0.05 





7.59±0.04 


7.59 ±0.04 


HIP 

Teff 

log 9 


logfFcI 

N(Fel) 

logfFcII 

N(Fell) 

Teff 

log 9 


logfFcI 

N(Fel) 

logepell 

N(Fell) 

18946 

4661 

4.54 

0.75 

7.63±0.09 

21 

7.63±0.11 

6 

4691 

4.64 

0.70 

7.56±0.09 

17 

7.56±0.21 

5 

13976 

5023 

4.51 

0.85 

7.69±0.07 

22 

7.69±0.11 

8 

5013 

4.61 

0.80 

7.64±0.08 

22 

7.64±0.08 

5 

19098 

5178 

4.59 

0.85 

7.62±0.09 

22 

7.62±0.06 

9 

5138 

4.54 

0.80 

7.63±0.09 

24 

7.63±0.06 

7 

16529 

5237 

4.51 

0.80 

7.64±0.10 

22 

7.64±0.08 

10 

5207 

4.51 

0.85 

7.62±0.08 

22 

7.62±0.05 

6 

19934 

5361 

4.57 

0.90 

7.64±0.06 

23 

7.64±0.06 

9 

5341 

4.57 

0.85 

7.62±0.08 

27 

7.62±0.04 

6 

20130 

5531 

4.45 

0.95 

7.66±0.09 

28 

7.66±0.07 

10 

5511 

4.55 

0.90 

7.62±0.06 

26 

7.62±0.07 

7 

20146 

5563 

4.33 

1.00 

7.67±0.08 

26 

7.67±0.08 

10 

5553 

4.43 

0.95 

7.62±0.07 

28 

7.62±0.04 

6 

19781 

5641 

4.35 

0.85 

7.69±0.09 

30 

7.69±0.08 

10 

5621 

4.40 

0.90 

7.60±0.06 

27 

7.60±0.05 

7 

19793 

5831 

4.31 

1.00 

7.73±0.08 

29 

7.73±0.09 

10 

5781 

4.41 

0.95 

7.64±0.08 

28 

7.64±0.05 

7 

20899 

5916 

4.31 

0.90 

7.66±0.06 

26 

7.66±0.08 

10 

5886 

4.31 

0.95 

7.61±0.08 

27 

7.61±0.04 

8 

19148 

6021 

4.37 

1.00 

7.61±0.06 

28 

7.61±0.09 

10 

5961 

4.37 

0.85 

7.62±0.08 

28 

7.62±0.04 

7 

22422 

6074 

4.34 

1.05 

7.66±0.09 

32 

7.66±0.08 

10 

6004 

4.34 

1.00 

7.63±0.06 

30 

7.63±0.01 

6 

21112 

6161 

4.36 

1.10 

7.60±0.09 

34 

7.63±0.04 

8 

6161 

4.26 

1.05 

7.60±0.06 

29 

7.60±0.03 

7 

22566 

6251 

4.30 

1.10 

7.69±0.08 

27 

7.69±0.09 

9 

6211 

4.30 

1.10 

7.66±0.08 

26 

7.66±0.07 

6 

average 




7.66±0.04 


7.66±0.03 





7.62±0.02 


7.62±0.02 
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Table 4. Abundance changes (A[X/H]) in a giant (HIP 20205) and in a main-sequence (HIP 19148) star in response to variations of T^g, logg and 


ASPL LIST 

El. 

T^eff + SO K 

Teff-SO K 

log p-tO. 10 dex logp-0. lOdex 

^-1-0.10 km s ' 

^-O.lOkms-* 





HIP 20205 




Fei 

0.03 

-0.01 

0.01 


0.00 

-0.04 

0.06 

Fell 

-0.04 

0.03 

0.05 


-0.05 

-0.02 

0.01 





HIP 19148 




Fei 

0.04 

-0.04 

-0.01 


0.01 

-0.02 

0.06 

Fell 

-O.OI 

0.01 

0.03 


-0.03 

-0.02 

0.01 

MASH LIST 





HIP 20205 




Fei 

0.03 

-0.03 

0.00 


-0.01 

-0.05 

0.04 

Fell 

-0.04 

0.04 

0.05 


-0.05 

-0.04 

0.04 





HIP 19148 




Fei 

0.04 

-0.03 

-0.01 


0.01 

-0.02 

0.02 

Fell 

-0.01 

0.02 

0.04 


-0.03 

-0.02 

0.02 


in Fe i and Fe it follow the same direction for both, the giant and 
main-sequence stars. 

We obtain a good agreement between the metallicity of gi¬ 
ants and dwarfs with this method. Moreover, since we are able 
to retrieve stellar parameters with a good agreement between Fe i 
and Fe it abundances, we do not see evidence for significant LTE 
departures either among the giants or the dwarfs according to 
this method. In Table |3l we find the metallicities obtained using 
MASH list are slightly higher than those using the ASPL list. 
The systematic offset of about ~ 0.04-0.06 dex is likely related 
to the different selection of lines, but can also be considered part 
of the internal errors of the analysis, as an accuracy better than 
0.05 dex can hardly ever be achieved without a line-by-line dif¬ 
ferential analysis. We further discuss this behavior in Section|5] 


4.3.1. Probing microturbulence velocities with 3D models 

Aside from the different techniques one can adopt to determine 
the microturbulence velocity, it is important to recall that this pa¬ 
rameter arises from a limitation of the classical ID model atmo¬ 
spheres into fully describing all the velocity fields present in the 
stellar photosphere. As a consequence, the inclusion of an ex¬ 
tra velocity field is required to describe the broadening observed 
in the li nes placed in the partly saturated regime of the curve of 
growth (IStruve & Elvevl[l9^1van Paradiislll972h . In practice, it 
is essential to use this parameter to obtain the same abundance 
for lines with small and large EWs. Thus, the optimal value of 
^ is obtained by imposing the absence of a correlation between 
the abundances and the EWs of a set of lines. This analysis, 
of course, depends on a good statistics of weak and moderately 
strong lines. This is usually a challenge in the simultaneous anal¬ 
ysis of giant and dwarf stars. 

One alternative to overcome the limitations cited above is 
to use 3D atmospheric models to predict Three-dimensional 
models treat convection in a physically consistent way, without 
the need of defin ing free parameters like f. This approach has 
been adopted by ISteffen et akl (1201 31 l2009l) , and it is based on 
the comparison between lines computed with spectral synthesis 
using 3D and ID models. Given a sample of spectral lines, the 
EWs computed from the 3D model with a fixed abundance are 
taken to represent the observation. Eor each line, the ID abun¬ 


dance is obtained by matching the “observed” 3D EW with the 
synthetic line profiles derived from the ID models (using exactly 
the same atomic line parameters as in the 3D synthesis). The 
best value of ^ for use with ID models is taken to be that which 
eliminates the correlation between line strength and derived ID 
abundance. In fact, the best ^ depends somewhat on the choice 
of the line list. Ideally, the lines should be insensitive to tem¬ 
perature fluctuations and have similar properties as the lines to 
be used for the abundance determinations. High-excitation Ee i 
lines and Ee ii lines are an obvious choice. 

We used 3D hydrody namic models taken from the CIEIST 


grid (iLudwig et al 


2009lrl computed with the CO^BOLD cod^^ 
(lErevtag et al.l2012 ). The metallicity of all selected models is so¬ 
lar since this grid does not have models available for stars more 
metal-rich than the Sun yet. Indeed, the differences between 
temperature stratification of 3D and ID models are expected to 
be small at slightly super solar metallicity. At this regime, the 
line opacities act toward heating the optically thin layers and, as 
a result, the mean temperature of the layer is close to radiative 
equilibrium. We do not expect the small metallicity differences 
to have any noticeable effect on the derived microturbulence ve¬ 
locities. 

The required ID h ydrostatic reference mo dels were com¬ 
puted with LHD code (ICaffau & Ludwig||2007h using the same 
stellar parameters and opacity scheme as the 3D models. Table|5] 
shows the stellar parameters of the computed models used for 
two giants, four subgiants, and six dwarfs, which were selected 
to be representative of the Hyades HR diagram. Eor giants and 
dwarfs, the models were selected to cover as well as possible the 
range of temperatures and gravities of our sample (e.g., Table|2l. 
Eor the subgiants, the models cover the cool part of the tempera¬ 
ture range expected for subgiants in the Hyades, but have slightly 
higher log g 3.50) with respect to what is expected from the 
cluster theoretical isochrone (log g ~ 3.0). Models for subgiants 
with log g smaller than 3.50 were not available in the grid. 

Eor computing the ^ values, we used 30 Ee i plus 7 Ee ii lines 
from the ASPL list and applied a method similar to method 3a 


^ An extended grid with respect to 2009. 

® conservative COde for the Computation of COmpressible COnvec- 
tion in a BOx of L Dimensions, L=2,3. http://www. astro.uu. se/ 
~bf/co 5bold_main.html 
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of IStefFen et alJ (l2013h . The selected lines are identified in Ta- 
h\dA2] Here again, lines with EW > 120mA were rejected to 
make this analysis compatible with the clipping criteria that were 
adopted in both Ml and M2. Briefly, for each line, we com¬ 
puted the equivalent width using the 3D model, W^u- Then, we 
computed for the very same lines a 2D curve of growth from 
the adopted ID reference model, VTiD(Aloge, ^), where Aloge is 
the abundance difference with respect to the original abundance 
used in the 3D spectrum synthesis. This grid allows us to find, 
by interpolation for given ^ value, Aloge, for each line i from 
the condition W^d. - fTio, - In other words, this abundance cor¬ 
rection, Aloge,, is the difference of the abundance computed in 
the ID model by fitting the equivalent width of the 3D line, and 
the true abundance used in the 3D spectrum synthesis. We com¬ 
puted Aloge, for a grid of microturbulence values ranging from 0 
to 2 km s ' with intervals of 0.1 km s *. The only exception was 
the giant with Teff = 5000 K, where the grid was ranging from 0 
to 3 km s“* in steps of 0.15 km s“*. 

For each ^ value, we plot Aloge, as a function of W^d (see 
Fig.E] left side, for an example). This graph illustrates the classi¬ 
cal concept of defining the microturbulence. We then determine 
the slope of the linear regression from each of the plots described 
above. These slopes are plotted against the corresponding micro¬ 
turbulence values, as shown in Fig.|3](right side). The best value 
of ^ for each star is taken to be that where the slope is zero in 
this curve. By following the procedure described here, we es¬ 
timated the microturbulence values for each of the model stars 
with parameters as listed in Table |5] 

Table 5. Microturbulence values computed from 3D models in corre¬ 
spondence with ID models. 


Star 

TeffK 

log g dex 

<^> 

Giant 

4477 

2.5 

0.89 

Giant 

4968 

2.5 

1.40 

Subgiant 

4582 

3.2 

0.84 

Subgiant 

4923 

3.5 

0.90 

Subgiant 

5432 

3.5 

1.13 

Subgiant 

5884 

3.5 

1.27 

Dwarf 

4509 

4.5 

0.71 

Dwarf 

4982 

4.5 

0.85 

Dwarf 

5488 

4.5 

0.87 

Dwarf 

5865 

4.5 

0.95 

Dwarf 

6233 

4.5 

1.02 

Dwarf 

6456 

4.5 

1.10 


We used the data listed in Table |5] to establish an empirical 
relation of the microturbulence as a function of effective temper¬ 
ature and surface gravity. We tested different functional forms 
for the calibration and adopted that which gave the best statisti¬ 
cal response. We searched for a calibration that shows no signifi¬ 
cant trend in the residuals distribution, but we have also assessed 
the quality of the fit by inspecting its correlation coefficient (R), 
the standard deviation cr of the fit, and the p-values for each 
term of the calibration. The best fit was found for the following 
equation: 

^ (kms"') = 0.998 h- 3.16 x lO^^'X - 0.253 Y 

-2.86 X IO^^'AT-h 0.165 T^ (2) 


where X = Teg - 5500 [K] and Y = log g - 4.0. The rms scat¬ 
ter of the residuals of this relation is 0.05 km s'*. However, we 
assume as the total uncertainty of the calibration a more con¬ 
servative value of 0.10 km s ’ because of the small number of 
stars use d in the fit. Although other calibrations in the litera¬ 
ture fe.g.. lBruntt et alJl2012tlEdvardsson et'^ll993L Feltzing & 
Gustafsson ll998h have considered a larger number of stars, their 
estimates of microturbulence are based on spectroscopic analy¬ 
ses with different line lists, while Eq. (I2]i presents a relation that 
reflects predictions from 3D models with a single line list, albeit 
for a limited set of stellar parameters. 

The results for the microturbulence v alues on Table |5] ar e 
consistent with those presented in Table 3 of ISteffen etal.l(l2013h . 
although these authors adopted a different selection of lines (Fe i 
lines with a lower excitation potential greater than 2 eV only) and 
a slightly different method (their method 3b). In this case, the 
microturbulence is given by the valu e that minimizes the scatter 
of the abundance corrections Aloge,. I Steffen et al.l (1201 3h prefer 
method 3b over 3 a, arguing that 3 a is more susceptible to details 
of the spectral line sample. Nevertheless, we decided to follow 
3a because it is closer to the usual analysis applied to determine 
^ in the literature, and thus facilitates the comparison with values 
derived in a classical spectroscopic analysis. 

We remark that our values display a trend for increasing f 
toward higher temperatures and (perhaps) lower gravities. How¬ 
ever, the limited number of stars tested here is not sufficient to 
provide more than a rough idea of the microturbulence behavior 
across larger areas of t he HR diagram. Mo re details about these 
trends can be found in ISteffen et al.l (l2013h . 

4.4. Method M2: The classical spectroscopic analysis 

The second method used to constrain the atmospheric parame¬ 
ters is the classical spectroscopic analysis based on the ioniza¬ 
tion and excitation equilibria of the Fe i and Fe ii lines. We want 
to investigate whether a consistent metallicity scale between gi¬ 
ants and dwarfs can be recovered through this method. Effective 
temperatures are calculated, forcing the Fe i line abundances to 
be independent of the excitation potential, i.e., forcing the exci¬ 
tation equilibrium. The microturbulence velocity is determined, 
forcing the Fe i abundances to be independent of the EWs. Sur¬ 
face gravities are calculated, forcing the lines of Fe i and Fe ii 
to produce the same abundance, fulfilling the ionization equilib¬ 
rium. As a consequence, the metallicity ([Fe/H]) is obtained as 
a spin-off of this procedure. All these criteria must converge in a 
fully consistent way and the final solution should be independent 
of the initial input parameters and iteration path. 

For this method, we use the same models as in Section 1431 
i.e, ID-FTE plane-parallel models (ATFAS9-HODFNEW). We 
also used the package WIDTH9 for the abundances computation, 
under some modifications to facilitate the handling of the in¬ 
put/output files. Extra IDLQ routines were written to optimize the 
calculation through the spectroscopic requirements mentioned 
before. First, the program computes the atmospheric model cor¬ 
respondent to the values of Tgg, log g, [Fe/H], and ^ given as 
initial guess. The code optimizes one parameter at a time, al¬ 
ways checking if the optimization of the previous parameter is 
still valid. If the previous optimization is still satisfied the code 
moves forward to the next parameter. If not, the code returns to 
the previous parameter and recalculates the optimization. The 
final set of spectroscopic parameters is obtained once all param- 

’’ IDL (Interactive Data Language) is a registered trademark of ITT Vi¬ 
sual Information Solutions. 
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Fig. 3. Microturbulence determination for the model of the dwarf with Teff = 5000 K and log g = 4.5 dex. Left: Aloge, versus the 3D equivalent 
width for ^ = 0.90 km s“* (symbols). The dashed line represents the linear regression to the data points. Right: slope of the linear regression as a 
function of microturhulence. The optimal microturhulence is derived hy the condition that the slope is zero. 


eters are optimized consistently. An example of the convergence 
of M2 is shown in Fig. |4] for the star HIP 20205, and for both 
ASPL and MASH lists. 

The internal errors of the M2 set of spectroscopic stellar pa¬ 
rameters were obtained as follows: for the Teff, we changed the 
slope of the linear regression in the [Fe/H] versus x diagram by 
its own 1 cr error. The error in temperature is the difference be¬ 
tween this new temperature and the previous best value. Simi¬ 
larly, for we applied the same idea, changing the slope of the 
linear regression in the [Fe/H] versus \og(EWIX) diagram. The 
log g error is estimated by changing this parameter until we ob¬ 
tain a difference between the Fe i or Fe ii abundances that is equal 
to the larger of their dispersions. For the metallicity, we adopted 
the standard deviation of the Fei abundance distribution. Ta¬ 
ble |6] shows the stellar parameters derived using M2 for ASPL 
and MASH lists. 

The metallicities derived by this method agree between gi¬ 
ants and dwarfs, within the presented uncertainties. However, 
the internal metallicity uncertainties of M2 are relatively larger 
than for Ml, when compared within the same line list. On aver¬ 
age, the internal uncertainties of M2 are about ~ 0.10 dex when 
using the MASH list and ~ 0.06 with the ASPL list. The dif¬ 
ference in the metallicity scale obtained for both lists are also 
larger when using M2. These differences are more significant 
for the giants and for the hottest dwarf of the sample. This result 
reflects how the line lists were assembled. The MASH line list 
was selected using the Sun as a proxy and has a larger number 
of transitions. This reduces the number of well-isolated lines, 
free from neighboring features, in the spectra of the giants. On 
the other hand, the ASPL list was chosen to have only isolated 
and uncontaminated transitions in the spectra of giants, which 
consequently reduced the number of lines used in the analysis. 
We discuss the quality of the line lists again in Section |5] On 
average, the difference between the atmospheric parameters ob¬ 
tained with the MASH and ASPL lists for M2 is -1-68+92 K for 
Teff, +0.12+0.24 for log g, -0.13+0.21 for^and+0.09+0.08 for 
[Fe/H]. These values are not significantly different from typical 
errors of spectroscopic analyses and from t ypical comparisons 
between multiple ana lyses methods (see, e.g.. lHinkel et al.l2014l : 
ISmilianic et al.ll^l4l) . 


4.5. Solar abundances and differential analysis 

We also computed the solar metallicity using our solar proxy 
spectra. This is useful to understand the behavior of the metal¬ 
licity scale derived with methods Ml and M2. As mentioned in 
Section [S] for the comparison with the giant stars we used the 
solar reflected spectrum of Ganymede, hereafter Sun HARPS. 
For the comparison with the dwarf stars, we used the solar spec¬ 
trum observed with UVES, hereafter Sun UVES. We decided 
to have two solar proxies to avoid any inconsistencies that may 
arise from the use of two different spectrographs and, therefore, 
different spectral resolution, instrumental profiles, or possible 
scattered light influence that affects a particular spectrograph. 

Solar abundances were derived applying the same two meth¬ 
ods presented before and for both MASH and ASPL lists. Eor 
Ml, in particular, we fixed the solar atmospheric parameters to 
5777/4.44/0.90 instead of ranging them within its expected er¬ 
rors. Eor the Sun, we prefer to keep fixed these parameters since 
its errors are too small to produce a noticeable difference in the 
metallicity determination. A 2cr clipping of the lines was also 
applied. Eor M2, the solar parameters were computed exactly as 
describe in Section l4[4l 

Table |7] shows the abundances of Eei and Eeii for our solar 
proxies, according to Ml and M2, using the ASPL and MASH 
lists. The solar atmospheric parameters found through M2 are 
very similar to the canonical values adopted as fixed in ML The 
mean differences between Ml and M2 for the solar parameters 
are 41 K for Teff, with minimum and maximum values of 16 K 
and 52 K; 0.01 dex for log g, with minimum and maximum val¬ 
ues of 0.01 dex and 0.02 dex; and 0.09 km s ' for f, with mini¬ 
mum and maximum values of 0.01 km s“* and 0.18 km s“*. 

Einally, we used the solar values of Table|7]as a reference and 
derived the abundances of Ee i and Ee it with respect to the Sun 
for the Hyades stars. The results of this differential analysis are 
shown in Table|8l for Ml and M2 and for both ASPL and MASH 
lists. Also given are the average metallicities obtained for giants 
and dwarfs according to each methodology. The internal disper¬ 
sions presented in the table were obtained by the squared sum of 
the internal uncertainty relative to the Ee abundance of the star 
and the Sun. 
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Table 6. Stellar parameters and iron abundances for giants and dwarfs derived according to method M2. The left side of the table shows the results obtained with the MASH line list and the right 
side the results obtained with the ASPL line list. 






MASH 







ASPL 




HIP 

Teff 

log 9 


logfFcI 

N(Fel) 

logfFcII 

N(Fell) 

Tcff 

log 9 


logepel 

N(Fel) 

logepcii 

N(Fell) 

20205 

4914±109 

2.88±0.07 

1.34±0.16 

7.62±0.11 

25 

7.62±0.07 

10 

4875±23 

2.71 ±0.06 

1.43±0.06 

7.53±0.07 

21 

7.52±0.06 

7 

20455 

5010±162 

3.03±0.09 

1.15±0.21 

7.72±0.13 

24 

7.72±0.11 

14 

4816±61 

2.55±0.07 

1.35±0.05 

7.50±0.06 

19 

7.51±0.07 

7 

20889 

4955±229 

3.36±0.11 

1.08±0.23 

7.86±0.16 

25 

7.86±0.12 

12 

4833±38 

2.74±0.08 

1.41±0.06 

7.59±0.09 

19 

7.59±0.07 

6 

average 




7.73±0.12 


7.73±0.12 





7.54±0.05 


7.54±0.04 


HIP 

Teff 

log 9 


logfFcI 

N(Fel) 

logfFcII 

N(Fell) 

TbB 

log 9 


logepei 

N(Fel) 

logepcii 

N(Fell) 

18946 

4815±141 

4.56±0.17 

0.91±0.10 

7.50±0.08 

20 

7.49±0.25 

10 

4813±104 

4.78±0.27 

0.54±0.20 

7.57±0.11 

18 

7.57±0.34 

7 

13976 

4900±87 

4.31±0.09 

0.70±0.19 

7.60±0.07 

20 

7.59±0.12 

9 

4915±42 

4.44±0.08 

0.63±0.17 

7.58±0.06 

20 

7.58±0.08 

5 

19098 

5120±75 

4.54±0.06 

1.00±0.19 

7.53±0.07 

21 

7.53±0.06 

9 

5025±62 

4.35±0.06 

0.93±0.13 

7.51±0.08 

25 

7.53±0.07 

7 

16529 

5090±86 

4.22±0.18 

0.83±0.13 

7.55±0.10 

23 

7.55±0.11 

11 

5100±48 

4.30±0.06 

0.83±0.10 

7.53±0.05 

20 

7.53±0.05 

6 

19934 

5330±59 

4.57±0.07 

0.81±0.14 

7.62±0.05 

21 

7.60±0.08 

10 

5200±45 

4.35±0.08 

0.81±0.11 

7.52±0.06 

25 

7.52±0.08 

7 

20130 

5590±104 

4.73±0.12 

0.98±0.17 

7.62±0.02 

17 

7.61±0.04 

8 

5505±40 

4.56±0.06 

1.01±0.07 

7.55±0.05 

25 

7.54±0.05 

6 

20146 

5500±92 

4.34±0.10 

0.58±0.12 

7.68±0.11 

27 

7.68±0.12 

11 

5540±40 

4.48±0.05 

0.96±0.08 

7.57±0.06 

27 

7.58±0.04 

6 

19781 

5695±53 

4.63±0.08 

0.78±0.12 

7.66±0.03 

17 

7.66±0.06 

8 

5625±41 

4.46±0.05 

0.96±0.07 

7.56±0.05 

27 

7.57±0.03 

6 

19793 

5790±88 

4.47±0.05 

0.80±0.20 

7.66±0.08 

28 

7.65±0.04 

7 

5710±34 

4.32±0.07 

1.13±0.07 

7.54±0.05 

26 

7.53±0.05 

7 

20899 

5855±105 

4.39±0.12 

0.85±0.17 

7.62±0.07 

28 

7.62±0.05 

8 

5885±47 

4.43±0.05 

1.09±0.07 

7.58±0.05 

26 

7.58±0.01 

5 

19148 

5985±91 

4.46±0.07 

0.75±0.15 

7.62±0.09 

30 

7.63±0.05 

9 

5970±50 

4.43±0.05 

1.09±0.07 

7.55±0.05 

26 

7.55±0.02 

6 

22422 

6029±104 

4.30±0.08 

0.57±0.28 

7.73±0.10 

32 

7.73±0.10 

11 

5975±53 

4.40±0.06 

0.94±0.09 

7.62±0.06 

30 

7.62±0.02 

6 

21112 

6135±113 

4.30±0.07 

0.83±0.24 

7.66±0.13 

36 

7.66±0.07 

11 

6095±44 

4.31±0.04 

1.16±0.09 

7.55±0.04 

24 

7.55±0.03 

7 

22566 

6343±152 

4.56±0.11 

1.24±0.14 

7.69±0.14 

33 

7.69±0.13 

11 

6010±59 

4.07±0.08 

1.16±0.09 

7.55±0.05 

24 

7.56±0.08 

6 

average 




7.62±0.06 


7.62±0.07 





7.S6±0.03 


7.56±0.03 
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Fig. 4. Fe I and Fe ii abundances as a function of the excitation potential (x), the logarithm of the reduced EW and the wavelength. The continuous 
red line is the regression between the quantities shown in each plot. The dashed line corresponds to the Feii final abundance derived according 
to M2. In all plots, the black asterisks correspond to the Fei lines and the red filled triangles correspond to the Feii lines. Left: spectroscopic 
analysis for HIP 20205 using ASPL list. Right: spectroscopic analysis for HIP 20205 using MASH list. 

Table 7. Solar abundances of Fei and Feii obtained with Ml and M2 for both MASH and ASPL lists. 


Ml 


MASH 



ASPL 





logfFel 

N(Fel) 

logeFdi 

N(Fell) 

logeFel 

N(Fel) 

logfFcII 

N(Fell) 

Sun UVES 

7.46±0.10 

36 

7.45±0.05 

14 

7.45±0.06 

31 

7.42±0.02 

5 

Sun HARPS 

7.45±0.07 

37 

7.44±0.06 

14 

7.45±0.05 

34 

7.43±0.02 

7 

M2 


MASH 



ASPL 





logfFel 

N(Fel) 

logfFdl 

N(Fell) 

logeFel 

N(Fel) 

logfFcII 

N(Fell) 

Sun UVES 

7.47±0.05 

22 

7.48±0.08 

15 

7.41±0.03 

25 

7.41±0.05 

6 

Sun HARPS 

7.50±0.05 

26 

7.49±0.07 

14 

7.42±0.05 

33 

7.42±0.02 

5 


5. Results and discussion 

5.1. The quality of the EWs 

One step in the spectral analysis that has a substantial influence 
on the final parameters and abundance is the measurement of 
the EWs. In this work, EWs were measured with the automatic 
code ARES. This cod e been used on many abundance analyses 
in the literature (e.g.. lAdibekvan et al.ll201^ iTabernero et al.l 
l2012h . Such kind of automatic codes are a fast and systematic 
way to compute EWs that minimize the subjectivity of manual 
measurements using, for example, IRAE. Nevertheless, it has 
been shown in the analysis of the Gaia-ESO Survey spectra that 
different groups measuring EWs in the same spectra, using the 
same code, can still find considerably differen t values (see the 
discussion in Section 6 of ISmilianic et al.f2014ll . This is the case 
mainly for two reasons: i) there are still some crucial free pa¬ 
rameters that need to be adjusted for the optimal measurement 
of EWs by automatic codes, in particular, for continuum fitting; 
and ii) quality control of the measured EWs is still important. 


as the codes just measure all possible EWs without recognizing 
potential local problems (e.g., unrecognized blends). Regarding 
the continuum placement, in general, there is no consensus on 
the best way to define the continuum, which can be done adopt¬ 
ing either a local or a global solution. 

It is not our aim in this paper to advocate in favor of choos¬ 
ing a global or local continuum normalization. We just stress that 
regardless of the type of normalization that was chosen, care is 
needed to evaluate whether the continuum solution was of good 
quality or whether it is affecting the measured EWs in a negative 
way. As ARES adopts a local continuum normalization, lines in 
crowded regions of the spectrum can have their EWs underesti¬ 
mated because of a low local continuum solution. This effect can 
become important in the spectra of giants, where many lines are 
stronger than in the spectra of dwarfs. Of course, if the line list 
is composed solely of well-isolated spectral features, the local 
and the global continuum solution should give similar results. 
However, this is not always the case, especially when analyzing 
stars in a wide range of T^g. Eor this reason, as stated before. 
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Table 8. Metallicity, with respect to the Sun, derived using Ml and M2. 




Ml 



M2 



MASH 

ASPL 

MASH 

ASPL 

HIP 

[Fel/H] 

[Fell/H] 

[Fel/H] 

[Fell/H] 

[Fel/H] 

[Fell/H] 

[Fel/H] 

[Fell/H] 

20205 

0.17±0.12 

0.18±0.11 

0.14±0.08 

0.16±0.06 

0.12±0.12 

0.13±0.10 

0.11 ±0.09 

0.10±0.06 

20455 

0.14±0.12 

0.18±0.11 

0.11 ±0.09 

0.13±0.06 

0.22±0.14 

0.23±0.13 

0.08±0.09 

0.09±0.07 

20889 

0.25±0.11 

0.26±0.12 

0.18±0.09 

0.20±0.06 

0.36±0.17 

0.37±0.14 

0.17±0.10 

0.17±0.07 

average - giants 

0.19±0.07 

0.21±0.05 

0.14±0.03 

0.16±0.03 

0.23±0.12 

0.24±0.12 

0.12±0.04 

0.12±0.04 

18946 

0.17±0.13 

0.18±0.15 

0.11±0.11 

0.14±0.21 

0.03±0.09 

0.01±0.26 

0.16±0.11 

0.16±0.34 

13976 

0.23±0.12 

0.18±0.12 

0.19±0.10 

0.22±0.08 

0.13±0.09 

0.11±0.14 

0.17±0.07 

0.17±0.09 

19098 

0.16±0.13 

0.17±0.08 

0.18±0.11 

0.21±0.06 

0.06±0.09 

0.05±0.10 

0.10±0.08 

0.12±0.09 

16529 

0.18±0.14 

0.19±0.09 

0.17±0.10 

0.20±0.05 

0.08±0.11 

0.07±0.14 

0.12±0.06 

0.12±0.07 

19934 

0.18±0.12 

0.19±0.08 

0.17±0.08 

0.17±0.08 

0.15±0.07 

0.12±0.11 

0.11 ±0.07 

0.11 ±0.09 

20130 

0.20±0.13 

0.21±0.09 

0.17±0.08 

0.17±0.07 

0.15±0.05 

0.13±0.09 

0.14±0.06 

0.13±0.07 

20146 

0.23±0.12 

0.24±0.09 

0.17±0.09 

0.17±0.04 

0.21±0.12 

0.20±0.14 

0.16±0.07 

0.17±0.06 

19781 

0.23±0.13 

0.24±0.09 

0.15±0.08 

0.18±0.05 

0.19±0.06 

0.18±0.10 

0.15±0.06 

0.16±0.06 

19793 

0.27±0.13 

0.28±0.10 

0.19±0.10 

0.22±0.05 

0.19±0.09 

0.17±0.09 

0.13±0.06 

0.12±0.07 

20899 

0.20±0.12 

0.21±0.09 

0.16±0.10 

0.19±0.04 

0.15±0.09 

0.14±0.09 

0.17±0.06 

0.17±0.05 

19148 

0.15±0.12 

0.16±0.10 

0.17±0.09 

0.20±0.04 

0.14±0.10 

0.15±0.09 

0.14±0.06 

0.14±0.05 

22422 

0.20±0.13 

0.21±0.09 

0.18±0.08 

0.21±0.02 

0.26±0.11 

0.25±0.13 

0.21±0.06 

0.21±0.09 

21112 

0.14±0.13 

0.18±0.06 

0.15±0.08 

0.18±0.04 

0.19±0.14 

0.18±0.11 

0.14±0.07 

0.14±0.05 

22566 

0.23±0.13 

0.24±0.10 

0.22±0.10 

0.25±0.07 

0.22±0.15 

0.21±0.15 

0.14±0.06 

0.15±0.09 

average - dwarfs 

0.20±0.04 

0.21±0.03 

0.17±0.02 

0.19±0.03 

0.15±0.06 

0.14±0.06 

0.15±0.03 

0.15±0.03 


we were careful to visually inspect and exclude measurements 
that could be underestimated because of continuum misplace¬ 
ment. This was particularly important in the MASH list, where 
the lines were basically chosen in the solar spectrum and, conse¬ 
quently, this list has spectral features that can be blended in the 
spectra of giants. 

We find that the EWs obtained with IRAF are systematically 
higher than the EWs obtained with ARES. We attribute this be¬ 
havior mainly to the different continuum normalization. The 
larger discrepancies were found for hot dwarfs and for the gi¬ 
ants. If we compare the EWs of all transitions, without exclud¬ 
ing those with suspicious normalization, the differences between 
IRAF and ARES reach up to 13.6 mA for the hot star HIP 22566 
and up to 7.4 mA for the giant HIP 20205. The use of all these 
measurements would certainly lead to an extra source of uncer¬ 
tainties. After the visual inspection and exclusion, the difference 
in EWs causes an effect on the abundances that is smaller than 
0.05 dex (see Section l4!2T i. This is similar to typical errors found 
in abundance analyses that are not line-by-line differential. 

We also noticed some difficulties of measuring EWs with 
ARES in hot stars where rotational broadening becomes to be 
important. After several tests with different ARES input options, 
we decided to limit our sample to those stars where we were 
confident that the EWs were well measured with this automatic 
tool. Therefore, we do not include stars hotter than ~6300 K. 

Finally, we also tested the influence of using spectra obtained 
with two different instruments on the EWs scale. To this end, we 
used a spectrum of the giant HIP 20455 observed with UVES, 
under the same configuration set as the dwarfs. We computed 
metallicities using MI for the two sets of EWs of this star. The 
mean differences were 0.03 dex for the ASPL list and 0.05 dex 
for the MASH list. Although the difference is slightly higher 
for the MASH list, it is still within the internal uncertainties of 
this method. Consequently, if there is any systematic error in 
the metallicity introduced by the use of different spectrographs, 


this is comparable to the internal errors of our analysis. There¬ 
fore, we concluded that the use of spectra from both UVES 
and HARPS instruments should not compromise the metallicity 
scale between giants and dwarfs. 

5.2. Comparison between the line lists 

The use of different line lists has already been reported as one of 
the possible sources of the discrepancies bet ween the metallicity 
scale of giants and dwarfs in open clusters (ISantos et al.ll2012L 
l2009h . However, in those works only the classical spectroscopic 
analysis was tested. Here, we are able to compared the spectro¬ 
scopic method (M2) with a set of parameters that is independent 
(method Ml). Thus, we can test the performance of the different 
line lists in these different approaches. 

When considering the analysis performed with the MASH 
list, the metallicity scale of giants and dwarfs agree better under 
M1. The difference between the average metallicity of giants and 
dwarfs is -0.02 dex with Ml and increases to 0.11 dex with M2. 
This increased difference seems to be driven by the giants. The 
uncertainties of the atmospheric parameters are larger for the gi¬ 
ants than for the dwarfs (Table [Q. Moreover, the metallicity and 
the surface gravity are substantially higher for two of the giants 
when M2 is applied with this list. On the other hand, with the 
ASPL list there is a better agreement between giants and dwarfs 
in both methods. The difference between the average metallicity 
of giants and dwarfs is -0.03 dex for Ml and -0.02 dex for M2. 

We attribute the increased scatter in the metallicity of the gi¬ 
ants seen when using the MASH list under M2 to the fact the 
lines were chosen from the solar spectrum. For this reason, they 
are likely more affected by blending features that were weak in 
the Sun, but are significant in the spectrum of a giant star. This 
result indicates that the use of fewer transitions, but that are care¬ 
fully chosen to avoid contaminating blends, is fundamental for 
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an accurate determination of atmospheric parameters and of the 
metallicity scale through the spectroscopic method. 

We also remark that the improved gf values that we adopted 
for the MASH list did not reduce significantly the scatter of the 
metallicities. Indeed, our values of metallicity dispersion are 
around 0.08-0.10 dex. T hese are similar to the ones reported in 
iMashonkina et al.l (1201 Ih . which are about ~ 0.09-0.11 dex. Per¬ 
haps the reduced number of lines in our analysis is not enough to 
reduce the internal statistical errors even improving the gf val¬ 
ues sources. For the list ASPL, we found an average dispersion 
on the metallicity of around 0.07 dex, which is smaller but still 
closer to the above mentioned values. 

We found a systematic difference of 0.04-0.06 dex between 
the abundances derived using MASH and ASPL lists, and higher 
for the giants when M2-I-MASH is applied. The metallicities de¬ 
rived with the MASH list are higher than the ones derived with 
ASPL list but also the total average temperature scale is hot¬ 
ter for the MASH list by about 20 K when using Ml and by 
about 70 K when using M2 (see Figure |5ll. These variations 
in Teff correspond to an increase in Fe abundances of approx. 
0.02-0.06 dex. In addition, the MASH list contains lines which 
are likely more affected by blending, thus we would expect that 
the cool giants show much larger abundances for the MASH list 
than for the ASPL list, which is the case. Thus, we attribute 
these systematic differences to the higher retrieved T^ff and to 
the presence of blended lines in the MASH list. More details 
about the comparison between Ml and M2 are presented in the 
next section. 

As mentioned at the beginning of this subsection, many 
works have demonstrated how abundances (and abundance dif¬ 
ferences or trends) depend on the adopted line list. From our 
analysis clearly emerges that a proper line list maximises the ro¬ 
bustness of the results. A similar conclusion has been reached by 
other authors who recently proposed to improve existing line list 
by customising it to the stars observed, using an empirical ap- 
proaclu_This_methodJias_been_aOTliedJo_cool_dwarfs_(Tsantaki 
et al. l2013h and to evolved stars (lAdibekvan et al.ll2015h . reach¬ 
ing substantial improvements with respect to what obtained with 
a not optimised line list. 

In terms of metallicity dispersion and compatibility of the 
metallicity scale between giants and dwarfs, the ASPL list shows 
a good performance for both methods tested in this work. This 
is likely related to the selection of well isolated line transitions 
which are equally appropriate for the analysis of giants and 
dwarfs. We, therefore, recommend the use of such list for this 
kind of analysis. 

5.3. Comparison between the methods 

In this Section, we compare the performance of methods Ml 
and M2 to investigate possible systematic effects between them. 
Figure |6] shows the difference between the atmospheric param¬ 
eters obtained with Ml and M2 as a function of the Ml param¬ 
eters. For each parameter, we performed a linear regression to 
check for significant trends. These trends would appear in case 
of considerable differences between the atmospheric parameters 
derived with each method. In case of a significant trend, the ra¬ 
tio between the slope and its own uncertainty should be less than 
0.5, i.e., X > 'IcT^. 

According to Fig.|6] only log g appears significantly differ¬ 
ent between Ml and M2, and only when using the MASH list 
(see the middle plot on the left panel of Fig.|6]l. On average, the 
spectroscopy surface gravities are larger when using M2 for both 
line lists, about 0.10 dex, but in the case of MASH the log g val¬ 


ues are significantly overestimated. This is probably caused by 
overestimated EWs, due to blending contaminants, which result 
in larger effective temperatures and surface gravities. The on 
the other hand, is higher on average for Ml than for M2. 

For the ASPL list, the plots shown in Fig. |6]do not indicate 
any significant difference between the parameters derived with 
Ml and M2. However, the Ml T^s values are, on average, ~ 47 K 
higher than the M2 ones. Although this is a small difference, it 
is probably responsible for the larger metallicities found with 
Ml in comparison with M2 (in absolute values, 0.05 dex for the 
giants and 0.06 dex for the dwarfs). 

The metallicities derived with both line lists show a better 
agreement when Ml is applied. This behavior can be more eas¬ 
ily noticed by inspecting Figure |5] The differences between the 
[Fe/H] scale derived with the two line lists with Ml are, on 
average, about ~0.05 dex, and constrained within the range of 
+0.1 dex (dotted lines) of this Figure. For M2, these differences 
are on average higher about ~0.10 dex, and present a larger scat¬ 
ter especially for the giants and for the hottest dwarf HIP 22566. 
Figure|5]shows that a good constraint in the atmospheric param¬ 
eters of a given set of star can produce a consistent metallicity 
scale for giants and dwarfs even when using line lists selected 
with different approaches, not optimized for the analysis of gi¬ 
ants. 

While the agreement between Ml and M2 is, on average, 
good, there are a few individual stars for which the differences 
are significant. The differences are about 100-150 K for the 
cooler stars (Tefr < 5300 K), with both line lists. Such discrep¬ 
ancies for the cooler stars do not compromise the overall agree¬ 
ment bet ween the two r.ff scale s, but this results are in agree¬ 
ment with lRamfrez et all (l2007h who argue that differences be¬ 
tween the spectroscopic and the photometric temperature scales 
are more significant for stars with < 5000 K. HIP 18946 and 
HIP 22566 are the stars that present the highest discrepancies 
between the effective temperatures of Ml and M2, reaching up 
to ~ 150 K for the MASH list. 

Finally, we tested the compatibility of the giants temperature 
scale obtained with the interferometric angular diameters and the 
IRFM. This last test was done in order too see if, in the absence 
of angular diameters for the giants, the IRFM scale would 
give similar results. As mentioned in Sect. 14.31 the agreement 
between the IRFM and the interferometric TeffS is excellent. The 
average absolute difference between them is -6+33 K. This re¬ 
sults indicates that, for similar giants, the IRFM could also 
be used as input values of temperature in ML 

In general, the metallicities derived using Ml are more con¬ 
sistent between giants and dwarfs for the two line lists tested 
here. This results seems to be a consequence of the well con¬ 
strained set of atmospheric parameters. The metallicities of gi¬ 
ants and dwarfs are only consistent when derived using M2 with 
the ASPL list. When using M2 associated with the MASH list 
we found systematic differences especially in the giants. Our re¬ 
sults confirm that, for the Hyades, the line list is a primary source 
of systematic differences on the metallicity scale of giants and 
dwarfs. 

5 . 4 . The Hyades metallicity scale 

In this Section we discuss the metallicity scale that we obtained 
for the Hyades. In general, the two methods tested favor the 
ASPL over the MASH list, as this last one seems to introduce 
systematic effects on the results, in particular for M2. Thus, we 
consider the results obtained with the ASPL list to be of higher 
quality. There is a good agreement between the metallicity of 
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Fig. 5. Comparison of the metallicity scale derived with lists MASH and ASPL as a function of the Teff obtained with ASPL list for Ml (upper 
panel) and M2 (lower panel). In each plot, the null, ±0.1 and ±0.2 abundance differences are represented by the continuous, the dotted and the 
dashed lines respectively. Red dots refer to the giants while black squares refer to the dwarfs. 
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Fig. 6. Comparison of the atmospheric parameters obtained with methods Ml and M2, as a function of the Ml parameters. The red lines are the 
linear regressions between the values. Left: parameters obtained with the ASPL line list. Right: parameters obtained with the MASH line list. 


giants and dwarfs when using this list for both methods. Ac¬ 
cording to Ml the difference between the average values of the 
iron abundances of giants and dwarfs is 0.03 dex. When M2 is 
considered, the difference between the metallicity of the giants 
and the dwarfs is 0.02 dex. These values are smaller than the in¬ 
ternal uncertainties of our methods and the typical errors found 
in classical abundance analysis (~ 0.05-0.10 dex). Therefore, we 


consider all the stars together, giants and dwarfs, analyzed with 
ASPL list to revisit the metallicity scale of the Hyades. 

We combined the average abundances of Feii for giants 
and dwarfs of each list to obtain the metallicity of the cluster. 
We choose Fe ii abundances because they are safeer abundance 
indicators than Fei lines, since Fei lines may suffer more of 
non-LTE effects. We found that the average metallicity of the 
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Hyades is 0.18+0.03 dex or 0.14+0.03 dex, according to Ml 
and M2, respectively, using the ASPL list. These values are 
compatible with each other within 2cr. For the MASH list, un¬ 
der the Ml analysis, we found a metallicity for the Hyades of 
0.21+0.04 dex. At this point, we reinforce two main aspects of 
our work. When investigating giants and dwarfs simultaneously 
using the classical spectroscopic analysis, it is recommended to 
use a line list especially suitable for the spectra of giants. In 
case the line list was selected based on the spectrum of a dwarf 
star, as the Sun, we recommend that it is more appropriate to 
use a method in the lines of Ml. This minimizes the systematic 
effects on the metallicity scale retrieved by the analysis. 


5.5. Comparison with other works 


We present here a comparison o f our results w i th sele cted pre- 
viousanalvses_ofll^ades_starsJPaulson_e^i2jl2003|^Schuler 
et al. (l2006l) and ICarrera & Pancinol (1201 Ih . These works were 
chosen because they have a larger number of stars in common 
with our analysis. As our final recommendation we select the 
results obtained w i th Ml using the ASPL list. 

Paulson et al. (l2003h analyzed a sample of 90 stars of the 
Hyades, of which 13 are present in our study. Those authors 
performed a differential spectroscopic analysis to derive the at¬ 
mospheric parameters of the stars. F or the sample in comm on, 
the average difference in Teff between fPaulson et al.l ( 2003 ) and 
this w ork is 58+59 K. The metallicity reported bv iPaulson et al.l 
(l2003h is 0.13+0.01 dex, which is within a 2cr range of the values 
found with Ml (0.18+0.03 dex), but it is in excellent agreement 
with our metallicity derived according the same spectroscopic 
method_^M2]_ofO. 14+0.03 dex. 

Schuler et al. (l2006l) analyzed the giants and dwarfs of the 
Hyades using a m ethod similar to Ml. The T^ff was ob tained 
with the IRFM bv iBlackwell & Lvnas-Gra^ (ll994fT9^ . Sur¬ 
face gravities and f were determined fitting synthetic spectra 
(ID+LTE) to the observed spectra. For the giants, the agreement 
between this work and our results is remarkably good for log g, 
^ and [Fe/H], and reasonable for Tf-g. The mean differences are 
76+12 K for T^g, 0.03+0.06 dex for log g, 0.06+0.05 km s“' 
for ^ and 0.02+0.01 dex for [Fe/H]. Among the eight dwarfs 
analyzed by this work, two were also included in our analysis. 
For these objects the agreement between the stellar parameters 
is good for T^^g and log g, but the differences increase up to 
0.40 km s ' for The overall metallicity for the Hyades recom¬ 
mended by these authors is 0.13 dex, which is also compatible 
within a 2cr range of our value retrieved with Ml and ASPL list 
(0.18+0.03 dex). 


Another analysis of the g iants of the Hyades was performed 
bv ICarrera & Pancinol (1201 ih . These authors determined the at¬ 
mospheric parameters via spectroscopy and photometry. The 
agreement betweenJhe_se^_of£arameter 2 _Ere 2 £ill£iiiLCarrera 
& Pancino (1201 ih and our work is reasonably good for T^g 
and good for [Fe/H], with average values of about -79+53 K 
and -0.03+0.03 dex, respectively. However, the surface grav¬ 
ities reported by these authors are compatible with those that 
we derived using M2, the classical spectroscopic analysis. The 
differences between the surface gravities is likely responsi¬ 
ble for the marginal agree ment between the m e tallici ty scales: 
0.11+0.01 dex derived bv ICarrera & Pancinol (1201 ih and our 
value of 0.18+0.03 dex. Although the value of metallicity re¬ 
ported by these authors is in excellent agreement with our value 
of 0.12+0.04 dex, which is the determination of the metallicity 
for the giants also derived using the spectroscopy method (M2 
and ASPL list). 


6. Conclusions 

We investigated the metallicity scale obtained in a simultane¬ 
ous analysis of dwarfs and giants. Inconsistencies between the 
metallicities of these types of stars have been reported before in 
the literature. Understanding whether these inconsistencies are a 
real or systematic effect of the analysis methods is important for 
advancing in a number of different areas of astrophysics. Of par¬ 
ticular interest is the comparison of metallicities between stars 
of different evolutionary phases hosting planets. 

As a test case, we chose a sample of giants and dwarfs in the 
Hyades open cluster. All stars in an open cluster are expected to 
share the same initial chemical composition. They are therefore 
optimal to test the consistency of analysis methods. 

We computed metallicities using two different methods and, 
for each method, using two different line lists. One line list was 
assembled based on the solar spectrum and with the aim of min¬ 
imizing non-LTE effects in the determination of the atmospheric 
parameters (MASH line list). The other was assembled with the 
specific purpose of analyzing giant stars, and therefore only in¬ 
cludes lines that are relatively free of blends in the spectra of 
these objects (ASPL line list). 

Analysis method Ml was based on atmospheric parameters, 
which are independent of spectroscopy. Effective temperatures 
were computed using interferometric angular diameters for the 
giants, and the IREM for the dwarfs. Surface gravities were 
computed with the help of theoretical evolutionary tracks. Eor 
use with Ml, we presented a new calibration of microturbulence 
based on 3D hydrodynamical models. In this way, one can also 
use an estimate of ^ that is independent of spectroscopic mea¬ 
surements. Method M2 was the classical spectroscopic method 
using Ee i and Ee ii lines to constrain the stellar parameters. 

We investigate the main steps in the analysis that may affect 
the metallicity scale of dwarfs and giants of the Hyades. Our 
careful evaluation of EWs has shown that differences in the con¬ 
tinuum placement can have an effect of about ~ 0.03 dex in the 
precision of the final metallicities. An effect of similar magni¬ 
tude (0.05 dex) can be introduced using spectra obtained with 
different instruments. Erom the stand point of the line lists, we 
found larger systematic differences between the metallicity scale 
for giants and dwarfs using MASH list, where the line selection 
criteria are not optimized for evolved stars. Additionally, we 
found a systematic difference of 0.04-0.06 dex between the two 
lists used in this work. Lastly, the difference on the metallicity 
scale for the giants and dwarfs using the two different methods 
(Ml and M2) is of about 0.04-0.06 dex, disregarding the results 
of the giants obtained with M2+MASH list. One interesting as¬ 
pect is that perhaps the use of the solar spectrum as a differential 
reference does not cancel out all these effects. Considering all 
these points, we conclude that the limit of precision of an abun¬ 
dance analysis that is not line-by-line strictly differential cannot 
be better than 0.03 dex. 

We show that with a careful determination of atmospheric 
parameters and a well -selected line list, it is possible to simul¬ 
taneously analyze giants and dwarfs and obtain consistency be¬ 
tween their metallicities. When attention to the line list is con¬ 
sidered, Ml and M2 produce results with a good agreement. The 
metallicity scale of the Hyades obtained with our preferred line 
list (ASPL) is 0.18+0.03 dex for Ml and 0.14+0.03 dex for M2, 
and these values are consistent among each other. At the very 
least, it is clear that assembling a line list well suited for the 
analysis of giants is mandatory to obtain consistency between 
the metallicity scales. When using the spectroscopy technique, 
the metallicity scale of giant and dwarfs using ASPL list agree 
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with previous spectroscopic analyses of the Hyades. However, 
Ml does show a more efficient capability to recover the metal¬ 
licity scale of giants and dwarfs regardless of the line list used. 
The results obtained with Ml for both MASH and ASPL lists are 
consistent among each other, indicating that a set of very well- 
constrained atmospheric parameters might compensate possible 
systematics from the features of the line list. Therefore, we be¬ 
lieve that Ml results are more robust. 

In the view of our results, we favor as our final metallicity 
for the stars in Hyades open cluster the value of 0.18+0.03 dex 
(based on Ml and the ASPL line list). For similar studies of gi¬ 
ants and dwarfs, in a more general approach, we suggest the use 
of Ml whenever possible. A simultaneous spectroscopic anal¬ 
ysis of giants an dwarfs, as done in M2, is only recommended 
under the use of a line list optimized for giant stars. 
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Table A.l. Atomic line data for the list MASH and iron abundances according to Ml. The data presented here corresponds to all the lines with 
EW < 120mA on the solar spectrum of Ganymede. Some of these lines were excluded as a result of a cr-clipping to produce the final abundances 
for the Sun. These are indicated with an * symbol in the last column of the table. 


Specie 

A (A) 

Alow (eV) 

logs/' 

EW (mA) 

logepc 

Clipped 

Fei 

4445.47 

0.09 

-5.410 

40.9 

7.470 



4574.72 

2.28 

-2.970 

54.9 

7.489 



4726.14 

3.00 

-3.250 

18.8 

7.588 



4808.15 

3.25 

-2.790 

26.9 

7.600 



4994.13 

0.91 

-3.002 

104.2 

7.275 



5197.94 

4.30 

-1.540 

36.9 

7.525 



5198.72 

2.22 

-2.113 

97.2 

7.396 



5216.27 

1.61 

-2.082 

117.7 

7.156 

* 


5217.40 

3.21 

-1.116 

111.3 

7.282 



5236.20 

4.19 

-1.497 

32.4 

7.298 



5247.05 

0.09 

-4.975 

66.5 

7.530 



5250.21 

0.12 

-4.918 

64.6 

7.462 



5285.13 

4.43 

-1.540 

27.9 

7.449 



5295.31 

4.42 

-1.590 

29.5 

7.518 



5379.58 

3.69 

-1.514 

61.5 

7.505 



5397.62 

3.63 

-2.528 

28.0 

7.702 

* 


5491.83 

4.19 

-2.188 

13.4 

7.445 



5517.06 

4.21 

-2.370 

17.5 

7.791 

* 


5522.45 

4.21 

-1.450 

42.9 

7.463 



5607.66 

4.15 

-3.437 

14.6 

8.699 

* 


5638.26 

4.22 

-0.770 

76.1 

7.394 

* 


5662.52 

4.18 

-0.573 

92.8 

7.402 



5679.02 

4.65 

-0.820 

59.0 

7.501 



5778.45 

2.59 

-3.430 

22.1 

7.409 



5807.78 

3.29 

-3.410 

8.5 

7.563 



5852.22 

4.55 

-1.230 

39.3 

7.470 



5855.08 

4.61 

-1.478 

23.8 

7.441 



5858.78 

4.22 

-2.160 

12.2 

7.383 



5916.25 

2.45 

-2.914 

54.1 

7.471 



5930.18 

4.65 

-0.230 

86.8 

7.335 



6065.48 

2.61 

-1.470 

116.7 

7.285 



6082.71 

2.22 

-3.571 

33.7 

7.456 



6105.13 

4.55 

-2.532 

12.8 

8.084 

* 


6151.62 

2.18 

-3.312 

50.0 

7.494 



6200.31 

2.61 

-2.405 

73.4 

7.497 



6213.43 

2.22 

-2.481 

82.6 

7.386 



6229.23 

2.85 

-2.805 

37.6 

7.379 

* 


6252.56 

2.40 

-1.727 

119.8 

7.437 



6421.35 

2.28 

-2.020 

103.7 

7.326 



6481.88 

2.28 

-2.985 

62.4 

7.497 



6518.37 

2.83 

-2.373 

59.0 

7.361 



6608.03 

2.28 

-3.930 

17.6 

7.428 


Fell 

4491.40 

2.84 

-2.710 

77.8 

7.468 



4508.29 

2.84 

-2.440 

84.8 

7.313 



4582.83 

2.83 

-3.180 

57.0 

7.455 



4620.52 

2.82 

-3.210 

50.3 

7.304 



5197.58 

3.22 

-2.220 

81.4 

7.315 



5264.81 

3.22 

-3.130 

44.3 

7.434 



5284.11 

2.88 

-3.195 

59.9 

7.525 



5414.07 

3.21 

-3.580 

25.6 

7.414 



5425.26 

3.20 

-3.220 

41.5 

7.429 



5991.38 

3.15 

-3.647 

29.5 

7.506 



6239.95 

3.89 

-3.573 

13.8 

7.635 

* 


6247.56 

3.89 

-2.435 

52.1 

7.512 



6369.46 

2.89 

-4.110 

18.6 

7.407 



6432.68 

2.89 

-3.570 

41.3 

7.442 



6456.38 

3.90 

-2.185 

61.5 

7.456 
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Table A.2. Atomic line data for the list ASPL and iron abundances according to Ml. The data presented here corresponds to all the lines with 
EW < 120mA on the solar spectrum of Ganymede. Lines indicated in the column 3D—^ list (f symbol) were also selected for the microturbulence 
calibration using 3D models. 


Specie 

d(A) 

Alow (eV) 

logs/' 

EW (mA) 

logepe 

3D-^ list 

Fei 

5044.21 

2.851 

-2.058 

73.9 

7.396 



5225.52 

0.110 

-4.789 

70.9 

7.482 

t 


5247.05 

0.087 

-4.946 

66.5 

7.501 

t 


5250.21 

0.121 

-4.938 

64.6 

7.481 



5651.47 

4.473 

-1.750 

18.6 

7.444 



5661.35 

4.284 

-1.756 

22.0 

7.369 

t 


5679.02 

4.652 

-0.750 

59.0 

7.431 

t 


5701.54 

2.559 

-2.216 

83.7 

7.504 

t 


5705.46 

4.301 

-1.355 

38.6 

7.373 

t 


5793.91 

4.220 

-1.619 

33.3 

7.435 

t 


5809.22 

3.883 

-1.710 

50.5 

7.516 

t 


5855.08 

4.608 

-1.478 

23.8 

7.441 

t 


5956.69 

0.859 

-4.605 

52.3 

7.532 

t 


6027.05 

4.076 

-1.090 

63.1 

7.397 

t 


6065.48 

2.609 

-1.530 

116.7 

7.345 

t 


6093.64 

4.607 

-1.300 

30.8 

7.395 

t 


6096.67 

3.984 

-1.810 

37.5 

7.455 

t 


6151.62 

2.176 

-3.299 

50.0 

7.481 

t 


6165.36 

4.143 

-1.460 

44.9 

7.454 

t 


6173.33 

2.223 

-2.880 

67.7 

7.484 

t 


6200.31 

2.609 

-2.437 

73.5 

7.532 

t 


6213.43 

2.223 

-2.520 

82.6 

7.425 



6240.65 

2.223 

-3.233 

48.5 

7.427 

t 


6252.56 

2.404 

-1.687 

119.8 

7.397 

t 


6265.13 

2.176 

-2.550 

84.9 

7.448 

t 


6270.23 

2.858 

-2.540 

52.2 

7.430 

t 


6430.85 

2.176 

-2.006 

110.0 

7.318 

t 


6498.94 

0.958 

-4.699 

45.8 

7.543 

t 


6703.57 

2.759 

-3.023 

37.7 

7.488 

t 


6705.10 

4.607 

-0.980 

46.5 

7.446 

t 


6713.75 

4.795 

-1.400 

21.1 

7.438 

t 


6726.67 

4.607 

-1.030 

46.7 

7.499 

t 


6750.15 

2.424 

-2.621 

74.3 

7.504 

t 


6810.26 

4.607 

-0.986 

48.9 

7.413 

t 

Fell 

5197.58 

3.231 

-2.220 

81.4 

7.415 

t 


5234.62 

3.221 

-2.180 

82.8 

7.388 

t 


5264.81 

3.230 

-3.130 

44.3 

7.461 

t 


5414.07 

3.221 

-3.580 

25.6 

7.414 

t 


5425.26 

3.200 

-3.220 

41.5 

7.453 

t 


6369.46 

2.891 

-4.110 

18.6 

7.407 

t 


6432.68 

2.891 

-3.570 

41.3 

7.442 

t 
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